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Abstract 

The static string-like solutions of the Abelian Higgs model coupled to dilaton gravity are 
analyzed and compared to the non-dilatonic case. Except for a special coupling between the Higgs 
Lagrangian and the dilaton, the solutions are flux tubes that generate a non-asymptotically flat 
geometry. Any point in parameter space corresponds to two branches of solutions with two different 
asymptotic behaviors. Unlike the non-dilatonic case, where one branch is always asymptotically 
conic, in the present case the asymptotic behavior changes continuously along each branch. 



PACS: 11.27.+d, 04.20. Jh, 98.80. Cq 



1 Introduction 

Of all the topological defects which may have been formed during phase transitions in the early 
tiniverse, cosmic strings ^ are those which have attracted most attention from a cosmological point 
of view. They were introduced into cosmology some 20 years ago by Kibble Q], Zel'dovich Q and 
Vilenkin Q , and were considered for a long time as possible sources for density perturbations and hence 
for structure formation in the universe. Indeed, the latest data from the BOOMERANG and MAX- 
IMA experiments ^, |j disagree with the predictions jlO| pT[ for the cosmic microwave background 
anisotropics based on topological defect models (see also ref. ]12|). This seems to point to the conclu- 
sion that if cosmic strings were formed in the early universe they could not have been responsible for 
structure formation. However, cosmic strings are still cosmologically relevant and enjoy wide interest in 
cosmology. 

The most common field-theoretical model, which is used in order to describe the generation of cosmic 
strings during a phase transition, is the Abelian Higgs model. This model is defined by the action: 




where TZ is the Ricci scalar, F^i, the Abelian field strength, $ is a complex scalar field with vacuum 
expectation value v and — — ieA^ is the usual gauge covariant derivative. We use units in which 
h = c = 1 and a "mostly minus" metric. 
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The simplest form of a cosmic string is just a direct generalization of the Nielsen-Olesen flux tube 
[p^ , which takes into account its gravitational effect. Because of the cylindrical symmetry, we use a line 
element of the form: 

ds^ = N'^{r)df - dr^ - L^{r)dip^ - K^{r)dz^ (1.2) 
and the usual Nielsen-Olesen ansatz for the +1 flux unit: 

$ = vf{r)e"^ , A^dxf" = ^(1 - P{r))dip (1.3) 

Since the flux tube is very concentrated along the symmetry axis, it is quite natural to expect that it 
will generate a spacetime geometry with the asymptotic behavior of the general static and cylindrically- 
symmetric vacuum solution of Einstein equations [ p^ , the so-called Kasner solution: 

ds^ = {krf^df - {krf'^dz^ - dr^ - (3^ {krf'^''-^\^ dip^ (1.4) 

where k sets the length scale while (3 represents the asymptotic structure, as will be discussed below. 
The parameters (a, 6, c) must satisfy the Kasner conditions: 

a + b + c = a^ + b^ + c^ = l (1.5) 

More information about the solutions can be obtained by inspection of the full system of field 
equations with the appropriate energy- momentum tensor, which we will not write here (see however 
Section 4 or refs.§, |l|l). Here we will only give a condensed summary. For cylindrical symmetry, 
the components of the energy-momentum tensor of the flux tube T^f have the property of — . 
This means that the solution will have a symmetry under boosts along the string axis, i.e., a = c. The 



Kasner conditions (1.5) then leave only two options, which are indeed realized as solutions of the full 
Einstein-Higgs system p^ . 

The standard conic cosmic string solution |6l E5[ is characterized by an asymptotic behavior given 



by ( |L4D with: 

a = c = , 5=1 (1.6) 

which is evidently locally flat. In this case, the parameter (3 represents a conic angular deficit |l8| , 
which is also related to the mass distribution of the source. 

In addition to the cosmic string solutions, there exists a second possibility: 

a = c = 2/3 , 6 = -1/3 (1.7) 



which is the same behavior as that of the Melvin solution |19[]. Eq.(1.7) is therefore referred to as the 
Melvin branch. Note however that the magnetic field in this solution inherits the exponential decrease 
(with r) of the original Nielsen-Olesen flux tube. Therefore, it is much more concentrated than in the 
original Melvin solution where it decreases only according to a power law. 

The main difference between the two branches lies in the Tolman mass, which is zero for the conic 
cosmic string solutions, but non-zero for the Melvin solutions [|l6). Moreover, the central magnetic field 
is generally larger for the Melvin solutions. 

The string-like solutions of the Abelian Higgs model are the simplest and most studied ones (see 
e.g. |l^-||l^, ||2^-||2^). It is, however, probable that when high energy corrections are taken into 
account, gravity is not purely tensorial. A minimal modification suggested by string theory is the 
introduction of a scalar degree of freedom, the dilaton 4>, turning gravity into a scalar-tensor theory in 
the spirit of Jordan-Brans-Dicke (JBD) theory. Studies of cosmic string solutions in the framework of 
JBD theory and its extensions already exist in the literature |2^, |2^ , and the typical characteristics 
of the string-like solutions are found to be quite different from those of the gauge strings of pure tensorial 
gravity. The main difference is the absence of asymptotically conic solutions due to the long range effect 
of the massless JBD field. Unlike the Brans-Dicke field, which is postulated not to couple to matter, 
the dilaton may do so and may further change the situation. This is represented by the following action 
dl, pi, which is written in the so called "JBD/string frame ||i": 



(1.8) 
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The parameter a serves as a general coupling constant between the dilaton and matter fields. Newton's 
constant is added explicitly in order to keep track of the dimensionality of the various fields. This makes 
the action in the "Einstein/Pauli frame |^" contain the standard Einstein-Hilbert Lagrangian: 



g |(7^ + 2V^0V''<; 



(1.9) 



It is related to the former by the metric redefinition: 



(B) 



(1.10) 



where gji^} is the metric tensor used in (1 



and gl^'' 



is used in (1.9). From now on we will use the 
Einstein frame for a clearer physical picture. 

This system was studied by Gregory and Santos j28j. Their analysis was done to first order in the 
parameter 7 = SttGv^ and is therefore valid only for weak gravitational fields. They concentrated in 
the one branch which contains for a — —1 asymptotically conic solutions. However, since for a ^ —1 no 
asymptotically conic solutions exist, there is no natural criterium to prefer one branch over the other, 
even if one is "more fiat" than the other. Moreover, the weak-field approximation breaks down both 
near the core of the string and far away from the string thus calling for a further investigation of the 
system. 

In this paper, we expand the study of cosmic strings in dilaton gravity in two directions. 

1) We analyze both branches and do not limit ourselves to weak gravitational fields. 

2) We modify the dilato n co upling to matter fields in such a way that it will be uniform and the original 
Einstein-Higgs system (1.1) is recovered in a certain limit. 



2 Field Equations and Asymptotic Solutions 



The field equations for the Einstein-Higgs-dilaton action ( |l.9D are straightforward to obtain and we 
give them already for flux tubes with cylindrical symmetry. Their structure becomes more transparent 
if we use dimensionless quantities. As a length scale we use I/VAw^ (the "correlation length" in the 
superconductivity terminology). We therefore change to the dimensionless length coordinate x — V Xv^r 
and we use the metric component V Xv^L which we still denote by L. We also introduce the "Bogomolnyi 
parameter" a = e^/A in addition to 7 = SttGu^, which has been already defined above. In terms of 
these new quantities we get, for fixed a, a two parameter system of five coupled non-linear ordinary 
differential equations (the prime denotes d/dx): 



I 2{a+l)<j> 7vr2 r fi\i / p 

if jJ^J^ + f e2(''+2)0 (1 „ /2) _ ^2(a+l), f^]f^O (2.1) 



r L 



^g2(a+l)0 y2p ^ Q ^2.2) 



(LNN'y 1 



7 (1-/2)2^ (2.3) 



^ f p2a0 J_ , 2{a+l),i> P , }_ 2{a+2)<t, _ f2^2\ /r, .^ 



N^L ^ 2a. 



+ /") + «6"^'£z^ + ^e^^"^'^'' (1 - ff) (2.5) 
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The sixth function, the metric component K{x), turns out to be equal to N{x) due to the high symmetry. 
We have also to keep in mind the existence of the constraint which comes from the (rr) Einstein equation, 
and gets the following form: 

TV' / L' N' 
— 2— H 

N \ L iV, 

/ ff2 p/2 p2 j2 1 \ 

+ -^fe2(a+l)0 i__^g2a0_^_g2(a+l)0:^_^_^g2(a+2)0(^_^2)2j (2J) 

The field equations are supplemented by the following boundary conditions, that should be satisfied 
by the scalar and gauge fields: 



/(O) = , lim fix) = 1 

x—^oo 

P(0) = 1 , lim P{x) = (2.7) 

X — ^oo 

These are just the usual flux tube boundary conditions borrowed from the flat space version. Moreover, 
regularity of the geometry on the symmetry axis a; = will be guaranteed by the " initial conditions" : 

L(0) = , L'(0) = 1 

iV(0) = 1 , N'{0) = (2.8) 

The presence of the dilaton requires additional conditions, which we impose on the dilaton fleld on the 
axis: 

0(0) = , 0'(O) = O (2.9) 
The condition (p'j O) = follows directly from the equations of motion and the constraint ( |2.6| ), using 



also eqs.( 2.7)-(|2.8|). As for the condition (j){0) = 0, notice that the equations of motion (2.1)-(2.6) are 



invariant under the following transformation: 

4> -> 4> + (f>o 

X xe""^" 

7 ^ ^^-2{a+l)<Po (2.10) 

where (j)o is a constant. Thus, we can trivially get the solution for 0(0) = (po from the solution for 
0(0) = by simple rescaling. 

Since we are looking for string-like solutions to the system, we may easily get the asymptotic behavior 
of the metric components from the assumption that the right-hand sides of the Einstein equations vanish 
exponentially fast, as in the case of pure tensorial gravity. Consequently, the system reduces to a simple 
set of equations, which are easily integrated to power law solutions: 

N{x) - KX^ , L{x) - Px^ , e'^(^) - Sx'^ (2.11) 

The 3 parameters A, B and C are subjected to two conditions, which may be viewed as a generalization 
of the Kasner conditions, eq.(|1.5D: 

2A + B^l , A^ + 2AB^C'^ (2.12) 

These two conditions leave a one-parameter family of solutions, which may be visualized by the ellipse 
in the B — C (or equivalently A — C) plane (flg.|l|) : 

^{B-^)' + 3C' = l (2.13) 

The 2 branches of the Einstein-Higgs system (the cosmic string branch and the Melvin branch) are 
represented in this picture by the top (po) and bottom (jjq) points of the ellipse, respectively. These 
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points correspond to an asymptotically constant dilaton (C = 0), and as we will show later, only the 
first is realized as an actual solution for the special value of the coupling constant a — 

The parameters A, B and C are related to the matter distribution through the following three 
quantities: The Tolman mass (per unit length) M, or rather its dimensionless representative m — GM: 

m^^ r dx N^L ( e^"-^ - }.e2ia+2)^ _ j2^2\ ^ 1 j.^ (^LNN') (2.14) 

and two others: 

3/2 



w — — 



I dx N^L (^e^"^^ ^ + e2(-+i)* + ie2(a+2)0 _ ^2^2^ ^ 



-i(lim(iV2L')-l) (2.15) 

Z X — >-oo 

which is related to the string angular deficit, and: 

D = jJ^°°dx N^L (^^e2(«+i)* + + cie"'^^ + ^e2(-+2)^ (1 - ff^ 

= lim {N'^Lcj}') (2.16) 

X — ^00 

which may be interpreted as a dilaton charge. Note that w is manifestly positive definite and D is 
positive definite for a > 0. The Tolm an mass is also non-negative since it is proportional to the power 
A which is non-negative due to ( ^.12| ). We further notice that, as in the non-dilatonic case w may 
be expressed in terms of the mass parameter m and the magnetic field on the axis, which we represent 
by the dimensionless parameter B: 

w^-B-m, where B = lim ( ^^^] (2.17) 
The parameters A, B and C are easily found to be expressible in terms of to, D and w or as: 

2to ^_ 2TO+l-g ^_ D ,2^g. 

6TO+1-6' 6TO+1-6' 6m-|-l-6' ^' ' 



The first of the two relations ( 2.12 ) is satisfied identically; the second translates in terms of to, D and 
;Bto: 

^ 4to(3to, + 1 - 6). (2.19) 



3 String-like Solutions 

Regarding the complexity and non-linearity of the system, it is quite evident that a detailed structure of 
the solutions can be obtained only numerically. We have performed a numerical analysis of the system 
and found that the solutions are paired, as in the non-dilatonic case, and have the characteristics that 
will be described below. The solutions have been found numerically by first discretizing the radial 
coordinate and then applying a relaxation procedure to the set of non-linear and coupled algebraic 
equations, T{f) = 0, in the set of scalar functions evaluated at the grid points, /. If the initial guess fi 
is good enough then the iteration obtained by /i+i — fi + A/,;, where A/,; is found by solving the linear 
equation: 

VT{h)^h = -^(/») (3.1) 

(Gauss elimination with scaling and pivoting will do), will converge towards a solution for the discretized 
variables. We constructed a starting guess satisfying the boundary and initial conditions and still having 
several free parameters. Changing gradually these parameters, we eventually stumbled into the region 
of convergence and a solution was found. Having obtained one solution of this type, the rest could easily 
be generated by moving around in the parameter space in sufficiently small steps. Representative plots 
of the solutions are given in figs.^j, || and ^. 
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The simplest picture emerges in the case a = —1. The two branches of solutions for this paticular 
value of a are particularly interesting. For fixed a, both solutions exist up to some critical 7; see fig. 5. 
One is an ordinary cosmic string, as in Einstein gravity with the additional feature of a non-constant 
dilaton. This solution has A = 0, B = l, C = independent of a and 7. The metric is very similar to 
the usual gauge string, but there are small deviations, due to the presence of the dilaton. It becomes 
identical to the metric of the Einstein gravity gauge string only in the Bogomolnyi limit (a = 2), where 



eqs.(|2T|)-(|2J) reduce to: 

P' = L{f - 1) (3.2) 
Pf 

/' = (3.3) 

L'-lpil-f) = l-l (3.4) 

with = 1 and (/> — 0, i.e. first order differential equations for L, f and P. A typical solution of this 
kind is depicted in figj|(a). 

The second branch in the case a = — 1 has for all values of a and 7 below the critical curve (fig. 5), 
C = — ^ and B = 0. Accordingly, A = ^. This solution is quite different from its non-dilatonic analog, 
the main difference being the asymptotic spatial geometry which is now cylindrical. 

At the critical curve (fig. 5), the two types of solutions become identical. Asymptotically the corre- 
sponding spatial geometry is conic with deficit angle Sip = 2tt. The critical curve itself can be approxi- 
mated by a power-law: 

7 « cia"^ (3.5) 

where ci « 1.66 and C2 ~ 0.275. This is quite, although not exactly, similar to the result obtained in 
tensorial gravity . 

The a and 7 - independence of A, B and C for a = — 1 is however the only case of such behavior. 
The generic one is such that for a given value of a, the powers of the asymptotic form of the metric 
components and the dilaton change with a and 7. Put differently, we have now a situation where the 
values oi A, B and C depend upon the three parameters a, a and 7. Therefore, the universal asymptotic 
form of the metric tensor, which corresponds to the two extremal points on the ellipse, is now replaced 
by a motion along the ellipse whose specific details depend on the value of a. Additional values of a 
were studied in detail and here we give few further results for two more special representative points, so 
in all there are: 

a = -73, -1, (3.6) 

These values were selected for the following reasons: Concentrating on the Einstein-dilaton-Maxwell 
system, a = — 1 corresponds to certain compactified superstring theories |3l|| , a = corresponds to 
Jordan-Brans-Dicke theory and a — — \/3 corresponds to ordinary Kaluza-Klein theory. 

In both cases a = — -s/S, a — the same two-branch picture emerges. One branch always contains 
the upper point of the ellipse (po) at the limit 7^0, such that for small 7 the "almost asymptotically 
conic" solutions of Gregory and Santos (28j are obtained. For large 7 the deviations from asymptotically 
conic geometry are more pronounced. At the same time a gravitational Newtonian potential appears, 
so this kind of string exerts force on non-relativistic test particles far outside its core. 

The other branch exhibits also an a and 7 dependence of the powers A, B and C, but has its 
7^0 limit at a-dependent points. This a-dependence can be understood if we note that the 7 — > 
limit corresponds to two types of asymptotic behavior. One is the asymptotically conic cosmic string 
behavior and the other is that of the dilatonic Melvin universe [ p9| ^ , which in our coordinate system 
corresponds to an asymptotic behavior with the following powers: 

2 - 1 2a , , 

^=^-^' 3.7 
-I- 3 -I- 3 -|- 3 

The situation therefore is the following. First we fix a. One branch (the "upper" or "cosmic string" 
branch) has the point po, i.e., A = C = 0, B = 1 a.s a. starting point for 7^0, and as a and 7 
change, the powers B and C move along the ellipse while A changes accordingly (see p. 12 ). The other 



branch (the "lower" or "dilatonic Melvin" branch) starts, for 7^0, with the above values eq.(3.7) for 
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the powers A, B and C and the motion along the eUipse is done in the opposite direction so the two 
branches approach each other and "meet" at a point in between. In the meeting point, the solutions 
not only have the same asymptotic behavior, but they become actually identical. 

There are no asymptotically conic solutions except for a — —I. Actually, this property can be 
understood analytically. For this, we note that a useful relation is obtained by adding the equation for 



the dilaton, ( |2.5| ) to the equation for the metric component N, (2.3) and integrating from to x, using 
the boundary conditions (on the axis): 

N^L{(t)' + N'/N) = 

7(a + 1) dxN^L Qe2(°+i)-^ (^^^ + /'^^ + e'"*^ + ie2(a+2)0 _ ^2)2^ (3 g) 

where we denote by x the integration variable. In the case a — —\, this equation is easily integrated 
again and gives in this case: 

N = e"-^ (3.9) 

It follows that the dilaton charge is minus twice the Tolman mass, D — —2m. This imposes an additional 
condition on the powers of the asymptotic metric components for any a, 7 which selects two points on 
the ellipse. The first is: 

A = C = , B = l (3.10) 

which is evidently locally flat. 

It is also easy to prove the converse, i. e., t hat there are no asymptotically conic solutions for a 7^ —1. 
In order to do it, we start again with eq.( |3.8| ) and notice also that for asymptotically conic solutions we 
need ^ = and B = 1, which yield C — 0. Thus, N'/N and </>' vanish asymptotically faster than 1/x 
and the left-hand side of eq. ( |3.8[ ) vanishes asymptotically. But the right hand side of this equation is an 
integral of a sum of 4 positive-definite terms and it can vanish only if a = — 1. 

It also follows that in this asymptotically conic solutions D = m ~ Q, while the central magnetic 
field is given by: 

B=l-K'f3=l-^(^l-^£^ (3.11) 

where Sip is the deficit angle. 

There is as usual a second point on the ellipse namely: 

^ = -C=l/2 , B = (3.12) 

As indicated above, this represents the same asymptotic behavior as the dilatonic Melvin solution (with 
a = —1). This solution has non-vanishing dilaton charge and Tolman mass. The central magnetic field 
fulfills the equation: 

B^l + 2m (3.13) 



4 Uniform Dilaton Coupling 

The model analyzed so far has a dilaton which couples (in the Einstein/Pauli frame) with different 
strengths to various matter fields. As a consequence, we do not return to the pure tensorial (Einstein) 
gravity in the limit of vanishing dilaton coupling. The only case of identical solutions of the two systems 
is the cosmic string solutions with a — ~1 and a = 2. It seems therefore natural to consider a uniform 
coupling of the matter fields to the dilaton. This uniform coupling will be described by the following 
action: 

S = J d'^xy/Wl (e^"^ Qi?^1>*D''$ - ^($*<D - - If.^F'^^^ + j^in + 2V^^W^)^ (4.1) 

Th e co rresponding field equations for the cylindrically symmetric case will be a simplified version of 
(131)- (H): 

(£!!;i^;££Z,,..*(a-/v5)/ = o ,4.2, 
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L 



/2p = 



(4.3) 



{LNN'Y 



2act> 



f P'- 



(4.4) 



-7 e 



2a4 



( P'^ p^ p 



£2 



(4.5) 



with the constraint: 



N' / L' TV' 

— 2 h — 

N \ L N 



pp2 



7 e M ^ 



2 



./2 



2aL2 



+ 7(1-/')' 



pr. 



2 j-2 



2aL2 2L2 



(4.6) 



(4.7) 



The asymptotic behavio r of t h e met ric components and the dilaton field are obtained as before and 
the results are identical to ( 2.11 )-( 2.1S ). 

Analysis of the new system shows that the solutions do not change much except in the case a = — 1, 
where an asymptotically conic solution does not exist any more, see fig. 6. Actually, there are no 
asymptotically conic solutions in this system for any combination of parameters and in this respect the 
solutions are more similar to the string-like solutions in JBD theory. It is quite simple to show that 
the only case where asymptotically conic solutions exist, in the case of this uniform coupling, is a — 0. 
The way to do it is to obtain an equation analogous to cq.(3.S) and to see that it can be satisfied by 
asymptotically constant N and only for a — 0. This is just the Einstein-Higgs system with a massless 
real scalar. 



5 Conclusion 

We have analyzed in detail, using both numerical and analytic methods, the cosmic string solutions in 
dilaton-gravity. Contrary to previous works, we did not limit ourselves to one branch in the weak-field 
approximation. We found it most natural to work in the Einstein frame, but all our results are trivially 
transformed into the string frame. 

We have shown that the solutions generally come in pairs. That is, any point in parameter space 
gives rise to two solutions with completely different asymptotic behaviors for the dilaton and for the 
geometry. Generally, neither of the two branches describe asymptotically flat solutions. This should be 
contrasted with the non-dilatonic case, where one branch always describes asymptotically flat solutions. 
Only for a very special coupling to the dilaton, which is related to certain compactified superstring 
theories, do asymptotically flat solutions appear. 

Our aim in this paper was to consider cosmic strings in the simplest dilaton-gravity models 
parametrized by a phenomenological coupling constant a. In the special case of fundamental string 
theory, which is related to our models for a = — 1, there is an infinity of curvature corrections (con- 
trolled by the reciprocal string tension a') and an infinity of string loop corrections (controlled by the 
string coupling which is proportional to e*^). However, since our solutions for a — —1 are everywhere 
non-singular and with finite dilaton, we expect that they will be only slightly modified in the full (yet 
unknown) quantum string theory. 

The asymptotic behaviors of the solutions in the two branches seemed at first to depend in an 
extremely complicated way on the parameters of the theory. However, a universal picture emerged, as 
explained in Section 3. 
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These results were obtained by coupling the dilaton "uniformly" to matter in the string frame. For 
comparison, we re-analyzed the whole problem using a uniform coupling in the Einstein frame which 
may seem to be more natural. Among other things, the correct tensorial (Einstein) gravity solutions 
now come out in the limit of vanishing dilaton coupling. But otherwise the solutions in the two models 
are quite similar. 
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Cosmic string branch with a=2, a— 1 and y=0.4 
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(a)" 

Dilatonic Melvin branch with a=2, a=-1 and 7^0.4 
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Figure 2: The solutions for a = —1, a = 2, 7 = 0.4. (a) Cosmic string branch, (b) Dilatonic Melvin 
branch. 
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Cosmic string branch with a=2, a=-V3 and 7=0.4 
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(a)" 

Dilatonic lyielvin branch with a=2, a=-V3 and 7=0.4 
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Figure 3: The solutions for a = — v^j a = 2, 7 = 0.4. (a) Cosniic string branch, (b) Dilatonic Melvin 
branch. 
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Cosmic string branch with a=2, a=0 and y=0.4 
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(a)" 

Dilatonic Melvin branch with oc=2, a=0 and y=0.4 
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Figure 4: The solutions for a = 0, a = 2, 7 = 0.4. (a) Cosmic string branch, (b) Dilatonic Melvin 
branch. 
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Curve where 5^=2n for a=-1 




Figure 5: The curve where Stp = 27r for the a 
ci w 1.66 and C2 « 0.275. 



= — 1 case. The curve is fitted to a power-law with 
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Cosmic string branch for the uniform coupling with a=2, a=-1 and 1^0.4 




Dilatonic Melvin branch for the uniform coupling with oc=2, a=-1 and ■^0.4 
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Figure 6: The solutions for uniform coupling with a=— l,a = 2, 7 = 0.4. (a) Cosmic string branch, 
(b) Dilatonic Melvin branch. 
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